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Abstract 

We use the modified quark-meson coupling (MQMC) model to study the composition profile of 
neutron star matter and compare the results with those calculated by quantum hadrodynamics 
(QHD). Both MQMC and QHD model parameters are adjusted to produce exactly the same 
saturation properties so that we can investigate the model dependences of the matter composition 
at high densities. We consider the possibility of deep kaon optical potential and find that the 
composition of matter is very sensitive to the interaction strength of kaons with matter. The onset 
densities of the kaon condensation are studied in detail by varying the kaon optical potentials. We 
find that the MQMC model produces the kaon condensation at lower densities than QHD. The 
presence of kaon condensation changes drastically the population of octet baryons and leptons. 
Once the kaon condensation takes place, the population of kaons builds up very quickly, and kaons 
become the dominant component of the matter. We find that the w-meson plays an important role 
in increasing the kaon population and suppressing the hyperon population. 
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I. INTRODUCTION 



Observation of neutron star properties such as mass, size and temperature provides us 
with important clues to the understanding of the state of matter at extremely high densities. 
In the 1970's, the maximum mass of the neutron star was calculated with the NN potentials 
available at that time 

Rasa 

and mean field models . Most of the calculations done 
in the 1970's resulted in stiff equations of state, and thus the maximum mass of a neutron 
star was predicted to be larger than 2M , where M is the solar mass. Only the Reid 
soft core potential yielded a soft equation of state and consequently a small maximum 
mass of a neutron star, 1.6M Q Recent observations of the masses of binary pulsars 
[?]], which are candidates of neutron stars indicate that the maximum mass of neutron 
stars are roughly around 1.5M , substantially smaller than most of the predicted values in 
the 1970's. (However, very recent observations seem to suggest possible existence of more 
massive pulsars in the range (1.8 — 2.O)M j^, though further confirmation is needed.) 
On the other hand, exotic forms of matter, i.e., matter consisting of degrees of freedom other 
than the nucleons have been proposed already many years ago. Some of the proposed exotic 
states of matter include those with creation of hyperons 10j, Bose-Einstein condensation 



(pions [llj or kaons [l2j), strange matter [13], and quark deconfinement [14|, [15|, [l6j . These 
exotic states seem to reduce the maximum mass of a neutron star close to the observations 
17, Q, Q, implying that exotic degrees of freedom seem to be needed to reproduce the 



observed masses of neutron stars. 

In this work, we consider the strangeness degrees of freedom by including both hyperon 
creation and kaon condensation in the neutron star matter. (It is the anti-kaon that matters 
here, but we simply refer to both kaons and anti-kaons as kaons for brevity.) The masses 
and energies of the hyperons and kaons in medium are sensitive to their interactions with 
the surrounding matter. In the meson-exchange picture, meson-hyperon and meson-kaon 
coupling constants can fix the strength of these interactions. The meson-hyperon coupling 



constants may be determined from the binding energies of hyperons in hypernuclei. 
meson-kaon coupling constants have been studied by using the kaon-nucleon scattering 

tl 



'he 



2d, 



2 if ] and kaonic atom data [20j. Recently, the magnitudes o 



matter have attracted much attention. Some calculations 



re kaon-nucleus potential in 



20 
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231 ] show that the real 



part of the K -nucleus optical potential Uk- is shallow (£/#■- ~ —50 MeV), but some other 
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calculations suggest that Uk- can be as large as about —120 MeV [2l|, |24| or even close to 



-200 MeV 



25]. 



Akaishi and Yamazaki predicted possible existence of deeply bound kaonic nuclei 26| , in 
which Uk- at normal density po was estimated to be about —120 MeV. Then, experiments 
at KEK claimed the observation of tribaryon kaonic nuclei, S° 27| and S + 28), which seem 
to suggest that K~ may be even more deeply bound than the theoretical prediction [26] 
(The former claim [271 ]. however, was withdrawn by the experimental group [29j). FIN- 
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301 ] and a BNL experiment with l6 0(K ,n) reaction [31] 



also reported distinct peaks. More recently there was a theoretical work which considered 
large kaonic binding energies and calculated widths of kaonic nuclear bound states 32] . The 
identities of these experimental peaks need to be studied further experimentally and theo- 
retically. However, in this work, we consider the possibility of deep optical potential of kaons 
in nuclei and explore the consequences in the composition profile of neutron star matter. 

In this work, for the description of dense matter we employ the modified quark-meson 
coupling (MQMC) model 33j. Nucleons and hyperons in the baryon octet are treated as 
MIT bags. The bag constant Bb and phenomenological constant Zb for a baryon B are 
fixed to reproduce the free mass of each baryon B. Coupling constants between (u, d)- 
quarks in the bags and (a, u, p)-mesons are adjusted to give us the binding energy per a 
nucleon Ej,/A = 16 MeV and symmetry energy a sym = 32.5 MeV at the saturation density 
Po = 0.17 fm~ 3 . Since the interaction between the s-quark and mesons are not well known, 
we adopt the standard quark counting rule and assume the s-quark is decoupled from (a, 
uj, p)-mesons. To take into account the interactions between s-quarks, we introduce cr* (980) 
and 0(1020) mesons following Ref. 18| for the baryon and Ref. 34( for the kaon. We also 
assume the kaon as a point particle. This treatment allows us to use Uk- as an input to fix 
the coupling constant between the cr-meson and the kaon, g a K- In our model the real part 
of the kaon optical potential at p = p can be written as Uk- = —{.9aK&(Po) + 9luI<u(po)), 
where cx(po) and cj(po) are the values of the meson fields at p . Using the value of g u x 
given by the quark counting rule, we can fix g a K for each given value of Uk-- Once the 
parameters of the model are fixed, the composition profile of neutron star matter can be 
obtained from the /3-equilibrium and charge neutrality. We find that the composition of 
neutron star matter changes dramatically depending on the value of Uk- ■ 

To investigate the model dependence of the results we also employ quantum hadrody- 
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namics (QHD) model 35] for calculating the composition of matter. The parameters of 
the QHD model are calibrated to produce exactly the same saturation properties as in the 
MQMC model. Our calculations show that the onset densities of the kaon condensation and 
the compositions of matter at high densities are substantially model dependent. In Sect. II, 
we introduce model Lagrangians and fix the model parameters. The results are discussed in 
Sect. III. Conclusions and discussions follow in Sect. IV. 



II. THEORY 



In this section we first briefly sketch the MQMC and QHD models by presenting the 
model Lagrangians. The models are calibrated so as to be consistent with each other at 
the saturation density by fixing the coupling constants of both models to produce exactly 
the same saturation properties; the saturation density, the binding energy, the symmetry 
energy, the nucleon effective mass, and the compression modulus. We then show how the 
physical quantities that will determine the composition of the neutron star matter can be 
obtained self-consistently. 



A. Models 



The model Lagrangian comprises the terms for the octet baryons, exchange mesons, 
leptons and kaons; C to t = C-b + £m + £-1 + £-k- Octet baryon, exchange meson and lepton 
terms in the mean field approximation can be written as 



B 



5> 



27 • d - m* B (a, a*) - 7 ( g wB uo + 
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+ -9pBT z po3 



if, 



B 



(i) 

(2) 
(3) 

where B denotes the sum over all the octet baryons (p, n, A, S + , S°, H°, and I 
stands for the sum over the free electrons and muons (e~, fi~). a, u and p mesons mediate 
the interactions between the non-strange light quarks (u and d). o* and <fi mesons are 
introduced to take into account the interactions between s quarks. Cb is of the identical 
form for both the MQMC and QHD models, but differs in the definition of the effective 
baryon mass m* B as will be shown below. 



MQMC 

In the MQMC model, a baryon is a composite system with quarks in a spherical bag, 
and its mass is given in terms of bag parameters an d q uark eigenenergy. The effective mass 



of a baryon in matter m B (a, a*) can be written as 



18 



33 



36. 



37 



38] 



m. 



(4) 



The bag energy of a baryon is given by 

E B 



EQ q 
~R 



(5) 



where Bb and Z B are the bag constant and a phenomenological constant for the zero-point 
motion of a baryon B, respectively. Q q - 



'x q + (Rm*) 2 , where m*(= m q — g%cr — g q a *<J*) 
is the effective mass of a quark whose free mass is m q . We take m q = for q = u,d and 
150 MeV for q = s. (Other choices of m q=s values do not make differences in the 



results 



391].) x q is determined from the boundary condition on the bag surface r = R, 

jo(x q ) =P q ji(x q ), (6) 



where (3 q = 
on density 



Qn—Rm* 



33 



,, . In the MQMC model, the bag constant Bb is assumed to depend 



38|. In this work, we use the extended form in Ref. [3] to include the 



contribution from a* as 



B B (cr, cr*) = i^Boexp 



4 



m B 



\ q=u,d q=u,d / ) 



(7) 



q=u,d q= 

where is the free mass of the baryon B, assuming that the a meson couples to u and d 
quarks only and that the a* meson couples to the s quark only. 
QHD 

In the QHD model, a baryon is treated as a point particle, and thus its effective mass is 
simply written as 



m B = m B - g^Bcr - g a *B<J 



(8) 



In order to reproduce the same saturation properties as obtained in the MQMC model, 
self-interactions of the a-field |4( 



(9) 



are added to Eq. (j2J) so that 



CT B = Cm - U^. (10) 

As mentioned above, the baryon and the lepton Lagrangians for the QHD model take the 
form given by Eqs. ([1]) and ([3]). 
Kaon 



The effective Lagrangian for the kaon may be expressed as 41j 



C K = D*K*D"K - m* K K*K, (11) 

where = <9 M + ig^K^^ — W<f>K ( t ) ^ + i\9pK^ '• P^- In this work we treat the kaon as a point 
particle in both MQMC and QHD models, and its effective mass is given by 

m* K = m K - g aK o - g a , K o* . (12) 

The equation of motion for a kaon is given by 

[D^ + m* K 2 ]K(x) = 0. (13) 

In uniform infinite matter, the kaon field K{x) can be written as a plane wave. Substituting 
the plane wave solution into the equation of motion, we obtain the dispersion relation for 
the anti-kaon 

UK = m* K - QukUq + g^K(f>0 ~ 9pK^P03- (14) 



B. Model parameters 
MQMC 

In the MQMC model, MIT bag parameters Bbo and Zb are determined to reproduce the 
free mass of a baryon B, m* B \ p= o = m B with the minimization condition ^§r~\ R - Rq = at 
a free bag radius Ro, which we choose as Ro = 0.6 fm. The bag parameters Bbo and Zb for 
the octet baryons are listed in Table [H 

Three saturation conditions p , E b /A, and a sym could determine three quark- meson cou- 
pling constants g„' d , g^' d an d g U p d ■, assuming u and d quarks to be identical in the isodoublet. 
The MQMC model, however, introduces an additional constant g'® in Eq. ([7]). Thus we 
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B 


m B (MeV) 


B 1 ^ (MeV) 


z B 


N 


939.0 


188.1 


2.030 


A 


1115.6 


197.6 


1.926 


£+ 


1189.4 


202.7 


1.829 


S° 


1192.0 


202.9 


1.826 


£- 


1197.3 


203.3 


1.819 


S° 


1314.7 


207.6 


1.775 




1321.3 


208.0 


1.765 



TABLE I: Bag constants Bbo and phenomenological constants Z B for octet baryons to reproduce 
the free mass of each baryon. The bag radius is chosen as Rq = 0.6 fm for all octet baryons, and 
the bare masses of quarks are fixed as rn u u\ = MeV and m s = 150 MeV. 



9% 


9% 


9 a 


9% 


m* N /m N 


K (MeV) 


1.0 


2.71 


2.27 


7.88 


0.78 


285.5 



TABLE II: The coupling constants between (u, <i)-quarks and (cr, u, /^-mesons in the MQMC 



model to reproduce the binding energy E^/A = 16 MeV and symmetry energy a 



sym 



32.5 MeV 



at the saturation density 0.17 fm 3 . m* N /m.N and K are the ratio of the effective mass to the free 
mass of the nucleon and the compression modulus at the saturation density, respectively. 

fix g^ ,d = 1, and adjust the remaining three constants to meet the three conditions. The 
resulting coupling constants are given in Table HT1 together with the ratio of the effective mass 
of the nucleon m* N /m^ and the compression modulus K. m* N and K are within reasonable 
ranges: m* N = (0.7 ~ 0.8)77^ and K = (200 ~ 300) MeV. 

The coupling constants between s-quarks and mesons cannot be determined from the 
saturation properties. In principle, experimental data from hypernuciei and kaon-nucleus 
scattering could be used to determine the coupling constants between s-quarks and mesons 
(for example, see Ref. 42j). However, these coupling constants are not well known yet, and 
for simplicity we assume that the quark counting rule holds and that the s quark does not 
interact with u and d quarks. Then we have 



n s — n s — n s — n u ' d — n u,d — 
9a — 9uj — Hp — 9a* —y<f> — u - 



(15) 



9aN 


9uN 


9pN 


92 (fm x ) 


93 


8.06 


8.19 


7.88 


12.139 


48.414 



TABLE III: The meson-nucleon coupling constants and the coefficients of the cr-meson self inter- 
action terms used in the QHD model. They reproduce the same saturation properties as in the 
MQMC model; p = 0.17 fm" 3 , E b = 16A MeV, a sym = 32.5 MeV, m* N = 0.78m N and K = 285.5 
MeV. 



To fix the meson-baryon coupling constants in the model Lagrangian, we also use the quark 
counting rule 

g P N = g P s = g P E = g% g P A = 0, 
l 

9<t>A — g<t>T, — -g^s — g$, 



(16) 



and the SU(6) symmetry 



/B 
9 a* 



V2, 
3.83, 



(17) 



The quark-meson coupling constants g% and g q p given in Table [TTl and the relations in 
Eqs. (15)-(17) determine all the meson-baryon coupling of the MQMC model. 
QHD 

In the QHD model, g a ^ and are adjusted to yield p and E b , and g p jy is fitted to 
produce a sym . g>2 and g% in JJ^ UD of Eq. are fixed to reproduce the same m* N and K 
values as listed in Table [TTJ for the MQMC model. The coupling constants determined in 
this way are given in Table ITTT1 In the MQMC model, meson-baryon coupling constants are 
obtained from the quark-meson coupling constants. On the other hand, in QHD meson- 
nucleon coupling constants provide the starting point for the determination of remaining 
other meson-baryon coupling constants. Once meson-nucleon coupling constants are fixed 
from the saturation properties, meson-hyperon coupling constants can be obtained by the 
quark counting rule (as in Eq. ( flBl ) and the SU(6) symmetry (as in Eq. ( TT7I) ). The coupling 
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U K - (MeV) 


-80 


-100 


-120 


-140 


-160 


g aK (MQMC) 


1.25 


2.01 


2.75 


3.50 


4.25 


g aK (QHD) 


1.26 


2.04 


2.82 


3.61 


4.39 



TABLE IV: g a K determined for several U K - values in the MQMC and QHD models. 

constants between strange mesons and hyperons can be obtained by combining the quark 
counting rule and the SU(6) symmetry, e.g., g^A = V%g u N/3 and g a *\ = V^g^N /3. 
Kaon 

There are 5 kaon-meson coupling constants in our models; g a K, g^K, 9 p k, 9u*k and g^K- 
g^K and g p x can be fixed from the quark counting rule: g U K = 9^ and g p x = gZ for the 
MQMC model, and g^K — 9l<jn/3 and g p x = g P N fo r QHD. (Obviously g p x from the MQMC 
model is the same as that from QHD. g U K{= 2.71) from the MQMC model is essentially the 
same as g U Ki= 2.73) from QHD.) g a * K may be fixed from / (980) decay [43j, and g^x from 
the SU(6) relation \[2g^ K = g n7rp = 6.04 [44J. g a * K and g^x thus fixed are 2.65 and 4.27, 
respectively. The remaining coupling constant, g a K, can be related to the real part of the 
optical potential of a kaon at the saturation density through Uk- = —(g<rK& + 5W^o)- gaK 
values corresponding to several values of Uk- are listed in Table IIVI for both MQMC and 
QHD models. 

Thus, out of 5 kaon-meson coupling constants g P K, g a *K-, and g^x are the same for both 
models. g w K^ are essentially the same for both models. g a K values are also very similar in 
both models for all Uk- values as seen in Table IV. Therefore, all the 5 kaon-meson coupling 
constants are practically the same for both MQMC and QHD models. 

C. Other quantities relevant to neutron star matter 

To obtain the composition of neutron star matter, we need to determine 16 unknown 
variables at each matter density, which include 5 meson fields (a,u, p,a*,4> ), 8 octet baryon 
densities, 2 lepton densities and the kaon density px- Five meson fields can be determined 
from their equations of motion: 

<° + To"'™ = E^C B (a) 2 -^±l £" ^^ ffl + (18) 
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m l*°* = Y,9**BC B (a*) 2J ^ 1 J —^— k 2 dk + g^ KPK , (19) 



^^(2Jb + 1)4/(6tt 2 ) - 9ukPk, (20) 



B 



= J2 9 ^ 2Jb + 1 ) fc i/( 67r2 ) +9*kPk, (21) 



m ?P03 = 9 P Bh B {2J B + l)4/(6vr 2 ) - ^-px, (22) 

B 

where J B and are the spin and the isospin projection, respectively, and k B is the Fermi 
momentum of the baryon B. In Eq. ffl8|) ^[/^ HD term needs to be there only for QHD and 
is not to be included in the MQMC model. C B (a) and C B (a*) are determined from the 
relations g oB C B {(r) = -^a and g a * B C B (a*) = ■ For QHD, C B (a) = C B (a*) = I. 



For MQMC, the explicit forms of C B (a) and C B (a*) are given in Ref. 18]. 
Charge neutrality condition of neutron star matter is expressed as 

^2 VbPb ~ P R ~ Pe ~ = °' ( 23 ) 

B 

where q B is the charge of baryon B and p B is the number density of B. Using the charge 
neutrality and the baryon number conservation conditions, one can fix two quantities, e.g., 
the density of the neutron and the electron. With these two variables fixed, /3-equilibrium 
conditions of the baryons give us the following 7 relations for the chemical potentials of 
p, A, S + , £~, E°, H~ and S° as 

p n = Pa = Ps° = Ps°, 

Pn + Pe = PS- = PS-, 

Pn - Pe = Pp = PX+, (24) 



where the chemical potential of baryon B is given by fi B = \Jk 2 B + m B 2 (<j, a*) + g wB ujQ + 
^0 + 9pBhBPo3- The chemical potential of a non-interacting lepton I is simply written as 



Pi = ^kf + mf. The /3-equilibrium condition for leptons 

Pe = Pa (25) 
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FIG. 1: Compositions of neutron star matter calculated from the MQMC (left panels) and the 
QHD (right panels) models. 



determines the density of muons. At a density where the condition 



(26) 



is satisfied, kaons are condensed and the kaon density becomes non-zero. Solving 
the Eqs. ( JT8H261) self-consistently and simultaneously, one can determine the 16 variables 
uniquely. 
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III. RESULTS 



Fig. [T] shows the relative populations, the ratios of the densities of octet baryons, leptons 
and K~ to the total baryon density, in the neutron star matter as functions of p/po up to 
p = 10po- The left panels show the results from the MQMC model and the right from the 
QHD model for U K - = -120,-140, and -160 MeV. (Figures for U K - = -80 and -100 
MeV are not shown here since they are not too much different from the one for Uk- = —120 
MeV particularly for QHD.) Figures from both models show that the onset density of the 
kaon condensation p crit becomes lower as \Uk- \ increases. 

To see how p crit changes depending on Uk-, let us consider Eq. (|2"6"|) . which determines 
p crit . Fig. [2J displays uj k and p n — p p , which are, respectively, the left and the right hand 
sides of Eq. ( 1261) (computed without including kaon condensation just for producing this 
figure). The left panel is from the MQMC model, and the right panel from QHD. At a 
density where the curve of ujk intersects with that of p n — p p , kaon condensation sets in. 
Among the coupling constants and meson fields that determine ujk, only g a x depends on 
Uk-. The a- meson contributes to ujk attractively, as can be seen in Eq. (|14p . Thus ujk 
becomes smaller for a larger |{7k--|, as shown in Fig. [21 

Figure [H also shows that as \Uk-\ increases from 120 MeV to 140 MeV p cr i t changes 
drastically in both MQMC and QHD models, but as \Uk-\ increases further above 140 
MeV, p cr it changes only moderately. This can also be seen from Fig. [21 As |£/j<--| increases 
from 120 MeV to 140 MeV the intersection between the curves for p n — p p and ujk moves 
rapidly (particularly for QHD), whereas when \Uk-\ increases further above 140 MeV the 
intersection shifts only a little to lower densities. 

Another common feature of the two models is that regardless of p cr it, once the kaon is 
created, the density of K~ piles up very quickly and overwhelms the population of the 



hyperons and even the nucleons. This behavior was also obtained by other authors 



41 



13, 



45l | . The reason can be partly attributed to the to- meson. The u;-meson term in the 



energy of K~ (in Eq. (14)) has a negative sign and is thus attractive, but it is repulsive for 
octet baryons. Fig. [3] shows the tu-meson is a dominant meson at higher densities in both 
MQMC and QHD models. Thus the tu-meson enhances the population of K~ but suppresses 
baryons, and the kaon density increases rapidly. In addition, due to the competition between 
the negatively charged hyperons and K~ in the charge neutrality condition, the negative 
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FIG. 2: Kaon energies u K for U K - = -120 MeV (dashed), -140 MeV (dotted) and -160 MeV 
(dot-dashed) are compared with /j, n — fi p (solid). At the densities where u>k and \x n — fj, p intersect, 
kaons start to condense. The left panel is for the MQMC model, and the right for QHD. 




FIG. 3: Meson fields calculated from the MQMC (left) and the QHD (right) models as functions 
of the matter density for U^- = —140 MeV. 

hyperons are highly suppressed and in some cases not even created at all as soon as the kaon 
condensation sets in. Positively charged hyperons, on the other hand, receive the opposite 
effects from the kaon condensation, and S + is created at lower densities as \Uk- \ increases 
more. The proton density is also enhanced by large abundance of K~, which facilitates in 
turn the enhancement of S + population through the chemical equilibrium condition of the 
positively charged hyperons in Eq. ( 1241 . 

Let us now discuss different aspects of the two model calculations. First, Fig. [1] shows 
that p crit from the MQMC model is lower than that from QHD. For U K - = -120, -140, 
and -160 MeV, Pcrit values are 5.9po, 3.8po and 3.0po in the MQMC model, respectively, 
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FIG. 4: Effective masses of the nucleons, kaon, and hyperons are plotted for Uk- = — 140 MeV. 

while they are 9.8po, 4.3po and 3.3po in QHD. Secondly, the MQMC model predicts a larger 
population of the kaon than QHD for a given Ujc- value. Figure [2] shows that ujk calculated 
from the MQMC model decreases more rapidly with density than ujk from QHD for each 
Uk- value. The curves for p n — \i v are more or less the same for both models at p < 4p , but 
at p > 4p Pn ~ Hp decreases faster in QHD. Thus the intersection and kaon condensation 
occur at lower densities in the MQMC model. This behaviour of the intersection in Fig. El 
is well reflected in the kaon condensation onset density p cr it in Fig. [lj Fig. [3] shows that the 
cr-meson field calculated by the MQMC model is larger than that calculated by QHD. A 
larger a- field in the MQMC model makes m* K and consequently lok smaller. On the other 
hand, as seen in Fig. El p n — Pp from QHD decreases faster with density at higher densities 
than that from MQMC. Thus the intersection of the ujk curve with the curve for p n — p p 
occurs at lower densities with the MQMC model. Therefore, p C rit is smaller in the MQMC 
model. 

Another model dependency of the results can be seen from the population of kaons, which 
is larger in the MQMC model. The effective mass of a kaon as a point particle is determined 
by cr and a* mesons through the relation m* K = uik — Qoko ' — 9a* k&* and is plotted in Fig. HI 
Since the a fields are larger in the MQMC model (as shown in Fig. [3]), the effective mass and 
the energy of a kaon are smaller in the MQMC model than QHD. Thus, kaon condensation 
takes place more in the MQMC model. 

Figures [3] and H] also show that even though cx-meson field from the MQMC model is 
larger than that from QHD as the densities increase, the reduction of the effective mass of 
baryons is smaller (or similar) with the MQMC model. If one could parametrize the effective 
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mass of baryons from the MQMC model in the form of m* B = m# — gaB{p)& — 5 , o-*b(c*)c* 
where ^o-s(cr) and g^B^*) are functions of a and cr*, respectively, the results in Fig. HJmight 
imply that (^(cr) and g^B^*) are decreasing functions with respect to the density. The 
rate of decrease is rather high since the product (go-s(cr) <t)mqmc is smaller than (or similar 
to) (g a B °")qhd while the cr-field value from MQMC is much greater than that from QHD. 
Such a decrease of g a B{&) in the MQMC model may be regarded as partial restoration of 
the chiral symmetry at high densities. 

We have calibrated both the MQMC and QHD model parameters to the same saturation 
properties. However, we find that the neutron star matter composition profiles from the 
two models are quite different and that they show significant model dependence. QHD 
assumes the baryons as point particles, whereas the MQMC model treats the baryons as 
MIT bags. Thus, the major difference between the two models is in the definition of the 
effective mass of baryons, m* B . The equation of motion for the cr-meson field is also different 
accordingly. m* B in QHD is a simple linear function of the cr-field, and the factor (7b (cr) in 
Eq (fl8l) . is a constant. In the MQMC model, m* B is a non-linear function of cr-field, and 
thus Cb(ct) is highly non-linear. When these non-linear m* B and Cb(ct) are expanded in 
powers of the cr-field, an infinite number of cr-field terms would appear. (Cubic and quartic 
terms are explicitly taken into account in the QHD model as in Eq. (Q.) Higher order terms 
can be interpreted as higher order contributions such as self interactions of meson fields, 
which are believed to be more important at high densities. But at high densities it can 
be questioned whether the non-linear terms of the cr-meson in the MQMC model account 
for the higher order effects properly and consistently. For instance, it is generally known 
that as the baryons come closer to each other the interplay of heavy mesons becomes more 
important. At high enough densities, their self interaction contributions may need to be 
included on the same ground as for the cr-meson, but the present MQMC model truncates 
the heavy meson terms at the leading order. 

It seems worthwhile to discuss at this point two more aspects of our results. The first 
one is the EoS and the resulting mass-radius relation of the neutron star. The second point 
is the dependence of our results on the £ hyperon interaction in matter, which is not well 
known yet. 

Let us first consider the EoS and the maximum mass of the neutron star. As the kaon 
(K~) appears and condensates, the number of negatively- charged hyperons decreases to 
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FIG. 5: Comparison of the EoS with and without kaons in QHD model. The Gibbs condition is 
used to treat the mixed phase. 

satisfy the charge neutrality. The decrease in the number of baryons will result in the 
reduction of the pressure and lead to softnening the EoS. Fig. [5] shows such a softening of 
the EoS due to the kaon condensation. In calculating the EoS of a system consisting of 
multicomponent substances, as in the case with the kaon condensation, Gibbs condition has 
to be employed for the proper description of the mixed phase. The curves in Fig. [5] are the 
results obtained with the QHD and Gibbs condition. As kaons appear, the EoS becomes 
considerably soft and the effect becomes more pronounced with a stronger attraction, i.e., 
for a larger \Uk-\ value. For the MQMC model, however, applying the Gibbs conditions 
do not give us a converging solution. Solving the 16 highly nonlinear equations together 
with Gibbs conditions double the number of equations to be solved, and the convergence 
could not be reached. It is not clear to us whether the convergence problem is due to 
numerical problems or due to non-linearity which can cause bifurcation or chaos. Therefore, 



we have used Maxwell construction for the MQMC model. (Some literatures 46] show 
that Maxwell construction is a good approximation to the Gibbs condition, but in some 
other literature [41[ it was emphasized that Gibbs condition produces significantly different 
results from those of Maxwell construction. Below we show that in our case the neutron 
star mass itself does not change much whether we use Maxwell or Gibbs conditions for 
QHD. Thus our use of Maxwell construction for the MQMC model may be considered as 
an acceptable approximation.) We solve Tolman-Oppenheimer-Volkoff equation to calculate 
the maximum mass of a neutron star. The results are shown in Tab. [V] where the central 
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MQMC (Mx) 


QHD (Mx) 


QHD (Gb) 


U K - 


p c /po M/Mq R 


p c /po M/Mq R 


p c /po M/M & R 


-120 


6.2 1.61 11.8 


6.1 1.50 11.4 


6.1 1.50 11.4 


-140 


4.6 1.53 12.8 


5.0 1.46 12.1 


5.0 1.45 12.1 


-160 


4.6 1.45 13.1 


4.0 1.32 12.7 


4.3 1.19 12.3 



TABLE V: The maximum mass of a neutron star M, the corresponding radius R and the density 
at the center of the star p c are listed for U^- = —120, —140, and —160 MeV. Uj<- is in units of 
MeV, and R in km. "Mx" and "Gb" refer to Maxwell and Gibbs conditions, respectively. 
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FIG. 6: Comparison of the population from the MQMC (left) with that of the QHD (right). The 
calculations are done with Uy, = +30 MeV and U K - = —140 MeV. 

density, the maximum mass, and the corresponding radius are listed for both MQMC and 
QHD models. The maximum mass calculated with QHD model is roughly 10 % smaller than 
that with MQMC model. For both models the maximum mass becomes smaller for a larger 
\U K \. The maximum mass calculated with MQMC and \U K -\ = 160 MeV is compatible 
with observation, while the maximum mass calculated with QHD and \U K -\ = 160 MeV 
becomes too small to be compatible with the observed values. However, this fact may not 
necessarily rule out the possibility of \Uk- \ becoming as large as 160 MeV because there are 
other possible mechanisms which are not included. 

Now let us consider the second aspect mentioned above. In the calculations made so far, 
we have assumed quark-counting rule in determining the hyperon-meson coupling constants. 
Experiments on A-hypernuclei indicate that quark counting is a good approximation of the 
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realistic interaction of A hyperons in nuclei, which gives the value of A optical potential 
at saturation density in the range —40 ~ —30 MeV. On the other hand, there is large 



ambiguity in the £ hyperon interaction strength. Ref. 47j shows that £ hyperon feels 



repulsion rather than attraction in nuclear medium. There are also experimental indications 



that the £ hyperon interaction is repulsive 48j|. If the £ interaction is indeed repulsive, 



the population profile of neutron star matter can change significantly from what is shown 
earlier in this work since we have used the quark counting rule. The number of £~ is closely 
related to the onset of K~ condensation since they compete with each other for the charge 
neutrality condition. To see the effect of possible repulsive nature of £ interaction on the 
kaon condensation, we have repeated the calculations with repulsive £ interaction. We first 
fit the coupling constants g'^ in MQMC and g^ in QHD so that the £ optical potential 
value at the saturation density equals to +30 MeV, and fix the remaining meson-£ coupling 
constants with the quark counting rule. The resulting population profiles with the kaon 
optical potential Uk- = —140 MeV are shown in Fig. [6j Compared to the quark-counting 
results in Fig. HJ the onset of kaon condensation occurs at slightly lower densities. This minor 
change happens regardless of Uk- value. However, the main features of kaon condensation, 
i.e., its onset density, fast increase of population and its dominance at high densities are not 
much affected by the change of £ interaction in nuclear medium. 

IV. CONCLUSIONS AND DISCUSSIONS 

Using the modified quark-meson coupling model, we have obtained the composition pro- 
files of neutron star matter, focusing on the effects of the strange particles of hyperons 
and kaons. Motivated by recent theoretical predictions of deeply bound kaonic states 2(| 



and the subsequent claims of the observations of interesting peaks found in KEK 27, 28], 

n n 

FINUDA[30j, and BNL[31[ experiments, large kaon optical potential Uk- was considered. 
By varying the value of Uk-, we have investigated how the onset density of the kaon con- 
densation and the composition of the stellar matter change. Employing the QHD model 
parameters which satisfy exactly the same saturation conditions as the MQMC model, we 
have investigated the model dependence of the results. 

We observed two common features from the two model calculations. First, a larger 
\Uk-\ produces a smaller onset density of the kaon condensation. This behavior is easily 
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understood from the relation between Uk- and g a K together with the role of g a K to the 
energy of the kaon ujk- Secondly, the number of kaons rapidly increases, and the number of 
negatively charged hyperons is strongly suppressed. This is due to the fact that the w-meson 
gives rise to attraction to K~ whereas it couples to baryons repulsively. 

Model dependence was also observed. The kaon condensation takes place at lower densi- 
ties in the MQMC model. The number of kaons is always larger with the MQMC model for 
given Uk- values. Larger cx-meson fields in the MQMC model can explain these behaviors. 
The differences in the results from the two models become more prominent at larger densi- 
ties. Growing discrepancies at higher densities have its origin partly in the effective mass of 
baryons in each model, which greatly affects the self-consistency condition of the cr-meson. 
The factor Cb(&) in the self-consistency condition of the cx-meson is highly non- linear in the 
MQMC model, which can be interpreted as an infinite number of cx-meson self-interaction 
terms. These higher order terms may require more proper and consistent treatment at high 
densities. 

An important issue in the dense matter physics is the restoration of the chiral symmetry. 
According to Ref . 49J , not only the mass but also the pion decay constant and meson- nucleon 
coupling constants decrease at a similar ratio at around the nuclear saturation density. In 
Ref. |50| the idea of scaling behaviour is applied to the neutron star matter using MQMC 
and QHD models with only nucleon degrees of freedom, and it is shown that the equation 
of state becomes stiffer when scaling effects are considered. This implies that if we include a 
scaling behaviour in our present models, it may ignite the onset of exotic states earlier than 
the present results which do not include a scaling. 

In the kaon sector, the coupling constants of a kaon and exchange mesons are currently 
an important issue. We took various values of the optical potential of K~ as an input to 
fix g a K- Other kaon- meson coupling constants are fixed from naive quark counting. It is 
known, however, that K + potential is repulsive with the magnitude Uk+ ~ 10 MeV at the 
saturation density [5l| . If Uk+ as well as Uk- is used as an input, g a K and g W K can be 
determined uniquely. For instance, if we take Uk- = —120 MeV and Uk+ = 20 MeV, then 
we get g a K = 2.041 and g^K = 4.187. g a K becomes smaller than the value listed in Table 
IV, while g^K becomes nearly twice of g^K fixed from the quark counting. Both a and u 
mesons contribute to the K~ energy attractively, but since the u meson becomes a dominant 
component at higher densities, taking into account of U K + can produce appreciably different 
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results. It may be interesting to see the effects of Uk+ on the kaon condensation. 

In our calculations, we have assumed the kaon as a point particle in both quark and 
hadron models. Comparison of the two models, however, shows that whether we treat a 
hadron as a bag (MQMC) or a point particle (QHD) can produce a significant difference. 
Therefore, it is worthwhile to treat the kaon g and compare the corresponding result 



with that of a point-like kaon. In Ref. [45( a kaon is treated as a bag in the framework of 
the QMC model, but no work has been done yet with the MQMC model. 

We assume m* K as a linear function of cr-field, but some authors employ a non-linear form 



mm 



m K = m K- 9aKm K a. (27) 



If we expand this expression in powers of a/rriK, we obtain 

m* K ~ itlk 



1 " 7,9*K— + 0{a 2 /m 2 K ) 



The leading order term of the cr-field has a factor 1/2, which is not present in Eq. (12). Due 
to the factor 1/2, the rate of decrease in m* K with density would be reduced by a factor 2, and 
this would shift the kaon condensation onset density to higher densities. This dependence 
on the kaon Lagrangian may be worthwhile to be studied. 
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